Interaction-induced adiabatic cooling and antiferromagnetism of cold fermions in 

optical lattices 
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We propose an interaction-induced cooling mechanism for two-component cold fermions in an 
optical lattice. It is based on an increase of the "spin" entropy upon localisation, an analogue of the 
Pomeranchuk effect in liquid Helium 3. We discuss its application to the experimental realisation 
of the antiferromagnetic phase. We illustrate our arguments with Dynamical Mean-Field Theory 
calculations. 
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Cold atoms in optical lattices Q] offer a promising labo- 
ratory for the study of strongly correlated systems, bring- 
ing quantum optics to have bearing on key issues in con- 
densed matter physics. Pioneering experiments on the 
Mott insulator to supcrfluid transition |2| have demon- 
strated the possibility of probing quantum phase tran- 
sitions between different ground-states of these systems. 
Recently, great progress has been achieved on cold Fermi 
gases as well, resulting in the production of molecular 
condensates in trapped gases 0, El an d the first 
imagingof Fermi surfaces in a three-dimensional optical 
lattice |8| . Controllability is one of the most remarkable 
aspects of these systems, with the possibility of tuning 
both the tunneling amplitude between lattice sites (t) 
and the on-site interaction strength ({/), by varying the 
depth of the optical lattice, and by varying the inter- 
atomic scattering length thanks to Feschbach resonances. 

In this letter, we consider fermionic atoms with two 
hyperfine ("spin") states in an optical lattice. When the 
lattice is deep and the scattering length is small (see be- 
low for a precise condition), a one-band Hubbard model 
is realized. The main physical effect studied in this paper 
is the possibility of cooling down the system by increas- 
ing the interaction strength adiabatically. As described 
below, this is due to a higher degree of localization -and 
hence an increase in spin entropy- as U/t or the tem- 
perature is increased. This is a direct analogue of the 
Pomeranchuk effect in liquid Helium 3. This mechanism 
relies on interactions and should be distinguished from 
the adiabatic cooling for non-interacting atoms in the 
lattice discussed in U 0] • The second main goal of the 
present paper is to study how this effect can be used in or- 
der to reach the phase with antiferromagnetic long-range 
order. For deep lattices (large U/t), the Neel tempera- 
ture is expected to become very low, of the order of the 
magnetic superexchange Jaf = 4< 2 /£7. Naively, it would 
seem that this requires extreme cooling of the gas. Here, 
we point out that the appropriate concept is actually the 
entropy along the antiferromagnetic critical line, and that 
at large U/t this quantity tends to a finite constant which 
depends only on the specific lattice. Hence, cooling the 



gas down to a temperature corresponding to this finite 
entropy per atom, and then following equal-entropy tra- 
jectories, should be enough to reach the magnetic phase. 
These physical observations are substantiated by theoret- 
ical calculations using in particular dynamical mean-field 
theory (DMFT) [H [12J , an approach that has led to im- 
portant progress on strongly correlated fermion systems 
in recent years. 

We consider the one-band repulsive Hubbard model: 
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(1) 



where i, j are site indices on the lattice, and a ="f, J. is a 
"spin" index associated with the two hyperfine states. 
The conditions under which two-component fermionic 
atoms in an optical lattice actually realize such a single- 
band lattice model will be discussed later. On an unfrus- 
trated bipartite three-dimensional lattice (e.g. the cu- 
bic lattice), with hopping between nearest-neighbor sites 
tij = t, and for one particle per site on average (half- 
filling), the physics of this model is rather well under- 
stood (see e.g. 13]). For temperatures above the Neel 
critical temperature T/v, the system is a paramagnet with 
an increasing tendency to Mott localization as U/t is in- 
creased (the Mott gap becomes of order U at large U/t). 
For T < Tn, the antiferromagnetic phase (Fig. dis- 
plays two-sublattice spin ordering and a doubling of the 
unit-cell. At weak coupling (small U/t), this is a spin- 
density wave instability with a weak modulation of the 
sublattice magnetization. In this regime, T/v is exponen- 
tially small in t/U, as a simple Hartree mean- field theory 
suggests. At strong coupling (large U/t), the low-energy 
sector of the model is described by a Heisenberg exchange 
Hamiltonian Jaf J2<ij> &i ' &j with Jaf = 4i 2 /£7. In 
this Heisenberg limit, T/v = 9 Jaf , with 9 a numerical 
constant depending; on the lattice (9 — 0.957 for the cu- 
bic lattice 13]). These two regimes are connected by a 
smooth crossover (which is equivalent to the BEC-BCS 
crossover at half filling) . The Neel temperature displays a 
maximum at intermediate coupling, as a function of U/t. 
This is illustrated by Fig. ^ on which we display our 
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FIG. 1: (Color online) Phase diagram of the half-filled Hub- 
bard model on the cubic lattice: antiferromagnetic (AF) and 
paramagnetic (PM) phase. Transition temperature within 
DMFT approximation (plain curve, open circles) and QMC 
calculation of Ref. |T^| (dot-dashed curve, squares). Dashed 
lines: isentropic curves (s=0. 4, 0.7, 0.75, 0.8), computed within 
DMFT. Dotted line: quasiparticle coherence scale Tp(U). 
The DMFT results were obtained with QMC (for T N ) and 
the IPT approximationfl if (for the isentropics) . The tran- 
sition curves are interpolations, continued at high U/t using 
the analytical expressions for the Heisenberg regime. 



calculation of Tn vs. U /t, using the DMFT approxima- 
tion on the cubic lattice and the quantum Monte-Carlo 
(QMC) Hirsch-Fye algorithm. DMFT overestimates T N 
by about 50% in the intermediate coupling regime, in 
comparison to the direct QMC calculations of Ref. |l3j 
on the cubic lattice (also displayed on Fig.^J. 

We now discuss how the entropy varies as the effective 
strength of the on-site interaction U / t is changed in the 
paramagnetic phase. Since all properties depend on the 
ratios T/t and U/t, we can consider that the hopping 
is fixed and that T and U are varied, or alternatively 
that both the temperature and coupling are measured in 
units of t, the natural unit of kinetic energy. Denoting 
by / and s the free-energy and entropy per lattice site, 
respectively, one has: s = —df/dT and df/dU — d, with 
d the probability that a given site is doubly occupied: 
d = (n^riii). We thus obtain: 



ds 
dU 



dd 
df 



(2) 



This equation can be used to discuss qualitatively the 
shape of the isentropic curves Tj = Ti(U) in the (U,T) 
phase diagram, along which s(T,(i7), U) — const. Taking 
a derivative of this equation yields: 



c{Ti) au - T *W T = T > 



(3) 



in which c = Tds/dT is the specific heat per lattice site. 
Fortunately, the temperature-dependence of the prob- 



ability of double occupancy d(T) has been studied in 
previous work by one of the authors [13. Hij and oth- 
ers [lfj. It was observed that, when U/t is not too large, 
the double occupancy first decreases as temperature is 
increased from T = (indicating a higher degree of lo- 
calisation), and then turns around and grows again. This 
is shown on Fig. [5] using DMFT calculations. This ap- 
parently counter-intuitive behavior is a direct analogue 
of the Pomeranchuk effect in liquid Helium 3: since the 
(spin-) entropy is larger in a localised state than when 
the fermions form a Fermi-liquid (in which s tx T), it is 
favorable to increase the degree of localisation upon heat- 
ing. The minimum of d(T) essentially coincides with the 
quasiparticle coherence scale Tp(U) which is a rapidly 
decreasing function of U (Fig. Q). This phenomenon 
therefore applies only as long as Tp > Tjy, and hence 
when U/t is not too large. For large U/t, Mott local- 
isation dominates for all temperatures T < U and sup- 
presses this effect. Since dd/dT < for T < T£(U) while 
dd/dT > for T > T£(U), Eq.© implies that the isen- 
tropic curves of the half-filled Hubbard model (for not too 
high values of the entropy) must have a negative slope at 
weak to intermediate coupling, before turning around at 
stronger coupling. In order to substantiate this behavior, 
inferred on rather general grounds, we have performed 
DMFT calculations of the isentropic curves, with results 
displayed in Fig^ The entropy s(T) was calculated by 
integrating the internal energy per site e(T) according 
to: s(T) = In 4 + e{T)/T - /~ dT'e(T')/T' 2 , which fol- 
lows from the thermodynamic relation dre = Tdrs. The 
DMFT equations were solved using the "iterated pertur- 
bation theory" (IPT) approximation ^lj (using, for sim- 
plicity, a semicircular density of states), and the internal 
energy was calculated from the one-particle Green's func- 
tion. 

It is clear from the results of Fig. [T] that, starting from 
a low enough initial value of the entropy per site, adi- 
abatic cooling can be achieved by either increasing U/t 
starting from a small value, or decreasing U/t starting 
from a large value (the latter requires however to cool 
down the gas while the lattice is already present). We 
emphasize that this cooling mechanism is an interaction- 
driven phenomenon: indeed, as U ft is increased, it allows 
to lower the reduced temperature T/t, normalized to the 
natural scale for the Fermi energy in the presence of the 
lattice. Hence, cooling is not simply due to the tunnel- 
ing amplitude t becoming smaller as the lattice is turned 
on. At weak coupling and low temperature, the cooling 
mechanism can be related to the effective mass of quasi- 
particles (tx 1/T£) becoming heavier as U/t is increased, 
due to Mott localisation. Indeed, in this regime, the en- 
tropy is proportional to T/Tp(U). Hence, conserving the 
entropy while increasing U/t adiabatically from (U/t)i to 
(U/t) f will reduce the final temperature in comparison to 
the initial one T f according to: T f /T, = T£(U f )/T£(Ui). 
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FIG. 2: Double occupancy d — (n^riii) as a function of 
temperature, for several values of U/t, calculated within 
DMFT(IPT). The initial decrease is the Pomeranchuk effect 
responsible for adiabatic cooling. 



At this stage, let us briefly discuss the validity of the 
DMFT approach, extensively used in the present work. 
In this approach, the lattice model is mapped onto a 
single-site quantum problem coupled to a self-consistent 
effective medium. This is an approximation, which be- 
comes exact only in the limit of infinite lattice coor- 
dination [rH . As a local approach, it underestimates 
the precursor antiferromagnetic correlations above TV, 
which will in turn quench the entropy and ultimately play 
against the cooling mechanism very close to TV- How- 
ever, as long as the correlation length is not too large, 
a local approximation should be accurate. Indeed, the 
existence of a minimum in d(T) has been confirmed by 
the calculations of Ref. using a different method, for 
a three-dimensional lattice, suggesting that the cooling 
mechanism discussed here is a robust effect. 

The isentropic curves in Fig.^suggest that interaction- 
induced adiabatic cooling could be used in order to reach 
the magnetically ordered phase. To explore this idea in 
more details, we focus on the entropy along the Neel crit- 
ical boundary sn(U) = s(Tn(U),U). At weak-coupling 
(spin-density wave regime), sn{U) is expected to be ex- 
ponentially small. In contrast, in the opposite Heisen- 
berg regime of large U /t, sn will reach a finite value sh, 
which is the entropy of the quantum Heisenberg model 
at its critical point, sh is a pure number which depends 
only on the specific lattice of interest. Mean-field theory 
of the Heisenberg model yields sh = In 2, but quantum 
fluctuations will reduce this number. We have performed 
a Schwinger boson calculation of this quantity, along the 
lines of fl8L Hif , and found that this reduction is of the 
order of 50% on the cubic lattice. How does sn evolve 
from weak to strong coupling ? A rather general ar- 
gument suggests that it should go through a maximum 
Smax> sh- In order to see this, we use again J2J) and take 
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FIG. 3: Phase diagram as a function of entropy. The dis- 
played curve results from a DMFT-IPT calculation (in which 
case sh = In 2), but its shape is expected to be general (with 
sh reduced by quantum fluctuations). 

a derivative of sn — s(Tn(U), U)), which yields: 

ds N _ c(T N ) dT N dd 
dU ~ T N dU dT lT=TN ( ' 

If only the first term was present in the r.h.s of this 
equation, it would imply that sat is maximum exactly 
at the value of the coupling where TV is maximum (note 
that c(Tjv) is finite (a < 0) for the 3D- Heisenberg model 
jUj)- However, in view of the above properties of the 
double occupancy, the second term in the r.h.s has a 
similar variation than the first one: it starts positive, 
and then changes sign at an intermediate coupling when 
Tp{U) — Tn(U). These considerations suggest that 
sn(U) does reach a maximum value at intermediate cou- 
pling, in the same regime where TV reaches a maximum. 
Hence, sn(U) has the general form sketched on Fig. [3] 
This figure can be viewed as a phase diagram of the half- 
filled Hubbard model, in which entropy itself is used as 
a thermometer, a very natural representation when ad- 
dressing adiabatic cooling. Experimentally, one may first 
cool down the gas (in the absence of the optical lattice) 
down to a temperature where the entropy per particle is 
lower than sh (this corresponds to T/Tp < sh/k 2 for 
a trapped ideal gas). Then, by branching on the opti- 
cal lattice adiabatically, one could increase U/t until one 
particle per site is reached over most of the trap: this 
should allow to reach the antiferromagnetic phase. As- 
suming that the timescale for adiabaticity is simply set 
by the hopping, we observe that typically H/t ~ 1ms. 

Let us now discuss the conditions under which two- 
component fermions in an optical lattice are accurately 
described by the Hubbard hamiltonian QJ (see also 
Q>E1)- The many-body hamiltonian is written in second- 
quantized form using as single-particle basis functions the 
Wannier functions associated with the periodic poten- 
tial V op t(r) — Vq sm 2 (nxi / a) (the lattice spacing is 
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FIG. 4: Spin-density wave and Heisenberg regimes as a func- 
tion of the depth of the periodic potential Vo and the scatter- 
ing length a a . The crossover between these regimes is indi- 
cated by the dotted line (U/t = 10), where Tn/t is maximum 
(other contour lines are also indicated) . In the shaded region, 
the one-band Hubbard description is no longer valid. Above 
the dashed line (U/A > 0.1) other bands must be taken into 
account and the pseudopotential approximation fails. Above 
the dashed-dotted line, non-Hubbard interaction terms be- 
come sizeable (td/t > 0.1, see text). 

a = A/2, with A the wavelength of the laser). The inter- 
action terms are obtained as matrix elements of the low- 
energy effective potential Vm^fx-f^) = 47r ^ a * S 3 (fl-r%) 
where a s is the scattering length. In general, this re- 
sults in a multi-band model which, besides the on-site 
Hubbard interaction, involves also more complicated in- 
teraction terms such as nearest-neighbor interactions, or 
density- assisted hopping terms of the form tdc\cjrii, with 
i,j neighboring sites. By explicitly computing these 
terms, as well as the one-body part of the hamiltonian, 
we examined under which conditions (i) the reduction to 
a one-band model is valid and (ii) these non-Hubbard 
interactions are negligible. This determines a domain in 
the (Vo/En,a s /a) plane (with Er = H 2/ K 2 /2ma 2 the re- 
coil energy) , which is depicted on Fig. Condition (i) 
requires in particular that the on-site Hubbard repulsion 
is smaller than the gap A between the first and the sec- 
ond band: U <C A. At large values of Vq/Eh, it can be 
shown that this is also the condition for our use of the 
pseudopotential approximation to be valid: a s -C with 
lo the spatial extension of the Wannier function of the 
first band. We found that the stricter condition of type 
(ii) originates from density-assisted hopping terms which 
should obey td <C t. We also displayed on Fig. 0] some 
contour lines associated with a given value of U/t. The 
one associated with U/t w 10 can be taken as the ap- 
proximate separatrix between the spin density-wave and 
Heisenberg antiferromagnetic regions. T^/t is maximal 
along this line, and T/v < 0.015Efj in the allowed re- 
gion. Thus adiabatic cooling is important to reach the 



AF phase. Since Vq and a s are the two experimentally 
tunable parameters, Fig.0]aims at summarizing useful in- 
formation for such experimental investigations. The de- 
tection of the antiferromagnetic long-range order might 
be achieved by spin-selective Bragg spectroscopy in order 
to reveal the doubling of the unit cell. The two hyperfine 
states could be distinguished by their Zeeman splitting, 
or using polarised light. A different method, which has 
been recently proposed jU| and investigated experimen- 
tally |2^|, is to use quantum noise interferometry. 

To summarize, we propose in this article an 
interaction-induced cooling mechanism for two- 
component cold fermions in an optical lattice. One 
possible application of this mechanism is in reaching the 
phase with antiferromagnetic long-range order. 
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